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Abstract: Compton scattering in the bulk accretion flow of the accretion column in magnetic cata- 
clysmic variables (mCVs) can significantly shift photon energies in the X-ray emission lines resulting 
from accretion shocks. In particular, Compton recoil can potentially broaden the 6.7 and 6.97 keV Fe 
Kq emission lines produced in the post-shock region, and contaminate the fluorescent 6.4 keV neutral 
Fe Kq line reflected off the white dwarf surface. We present nonlinear Monte Carlo simulations demon- 
strating these effects, and we discuss the interpretation of observed Fe Kq linewidths in mCVs in light 
of these new results. The implications for other accreting compact objects are also discussed. 



Keywords: accretion: acretion columns 
- stars: white dwarfs - X-rays: spectra 



line: profiles - scattering: electron - stars: binaries: close 



1 Introduction 

Magnetic cataclysmic variables (mCVs) are low-mass 
close binaries with a Roche-lobe overfilling red dwarf 
transferring material to a magnetic white dwarf (see 
e.g. Warner 1995; Cropper 1990). The two main sub- 
classes of mCVs are the polars (AM Herculis systems) 
and the intermediate polars (IPs). In a polar, the 
white dwarf has a sufficiently strong magnetic field 
(B ~ 10 - 200 MG) such that the white dwarf and its 
companion star are locked into synchronous rotation. 
The magnetic field also prevents the formation of an 
accretion disk, causing the material from the red dwarf 
to be accreted onto the white dwarf directly along its 
magnetic poles. Although the magnetic field of a white 
dwarf in an IP is weaker, allowing an accretion disk to 
form, the inner part of the disk is truncated by the 
white-dwarf magnetic field. In both polars and IPs, 
the accretion fiows onto their white dwarfs are con- 
fined by the magnetic field, and accretion columns are 
formed, channelling the accreting matter to the white- 
dwarf magnetic polar regions. 

An accretion shock often develops in the accre- 
tion column such that the supersonic accretion flow 
becomes subsonic in the downstream, post-shock re- 
gion, before settling onto the white-dwarf surface. The 
shock heats the accreting matter to temperatures 



3 GM^finiH 
8 (iiw + Xs) 



(1) 



where Mw and J?w are the white-dwarf mass and ra- 
dius, fi is mean molecular weight, and Xs is the shock 
height, which is typically ^ _Rw (see, however. Crop- 
per et al. 1999 for cases with large shock heights). 
For white-dwarf masses of 0.5 — l.OM© and typical 
mCV parameters, the shock temperature is kTs ~ 10 — 
40keV (~ 1 - 4 X 10* K). The hot plasma in the post- 



shock region is cooled by emitting bremsstrahlung X- 
rays and optical/infra-red cyclotron radiation (Lamb 
& Masters 1979; Wu 2000 and references therein). 

At a plasma temperature of kT ~ 1 — 50 keV, 
heavy elements such as Iron (Fe) are highly ionised, 
resulting in an abundant population of H- and He-like 
ions (see Fujimoto 1998 and Wu, Cropper & Ramsay 
2001 for the ionisation structure of post-shock flow in 
mCVs). The 6.97 keV Lyman-a line is emitted due to 
the 2p Is '^S transition in the H-like Fe XXVI 

ions. The 6.70 keV He4 line is emitted from the He- 
like Fe XXV ions for the 2p ^Pi Is ^^o transi- 
tion; and the 6.67 and the 6.68 keV He5 lines, for the 



2p ^Pi 



Is ^5*0 and the 2p ^P2 Is ^So transi- 



tions respectively. The 6.64 keV He6 line, due to the 
2p ^Si Is ^So forbidden transition, is generally sup- 
pressed in the shock-heated regions of mCVs, because 
of a high (> 10^'' cm~'^) electron number density. 

The natural widths of these Fe Kq lines are small. 
Realistically, the lines will be Doppler broadened by 
the bulk motion of the emitters in the flow, and also by 
scattering processes. X-ray observations of mCVs have 
indeed rev ealed substantially broadened Fe Kq emis- 
sion lines (iHellier et al.lll998l : iHel licr & Mukai 2004), 
though flow-velocity induced Doppler shifts are yet to 
be verified. The line photons may undergo a large 
number of resonant scatterings with the ions, but this 
requires the scattering region to have small velocity 
gradients (See Terada et al. 2001 for a detail study 
of resonant scattering in mCVs). Electron (Comp- 
ton) scattering can shift line photon frequencies signif- 
icantly and thereby allow the photons to escape res- 
onance trapping. It also can cause substantial line 
broadening. 

A previous study (Wu 1999) proposed that cold 
electrons in the upstream pre-shock bulk flow of mCVs 
can scatter a significant number of line photons emit- 
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ted from the shock-heated region, and any 6.4 keV flu- 
orescent Fe hne from the white-dwarf surface could be 
contaminated by the down-scattered 6.7 keV He-hke Fe 
hne photons. A recent study (Matt 2004) showed that 
such scattering could be important in that it modifies 
the X-ray polarization properties. In this paper, we in- 
vestigate the effects of electron (Compton) scattering 
in the accretion column of mCVs by means of numer- 
ical simulations. We focus on its effects on the H-like, 
the He-like, and the neutral Fe Kq lines. We employ 
a non-linear Monte Carlo algorithm (CuU en 2001a,b). 
which improves on other similar codes ([Stern et alJ 
Il995 : Hua 199711 by taking into account nonuniformity 
in the velocity, density and temperature profiles in the 
scattering region. The paper is organised as follows: 
in §2, we briefly discuss various line broadening mech- 
anisms; in §3 we outline the geometry and physical pa- 
rameters of our model for the line-emitting and Comp- 
ton scattering regions in mCVs, and we summarise the 
Monte Carlo method used in our simulations; results 
are presented and discussed in §4; and conclusions are 
given in §5. 



m 




2 Broadening of Fe Ka Lines 
from mCVs 

The post-shock region of mCVs are mostly optically 
thin to bremsstrahlung X-rays, and the coronal ap- 
proximation is usually applicable for line emissivity 
calculations (e.g. Fujimoto 1998; Wu, Cropper & Ram- 
say 2001). The thermal Doppler widths of the Fe Ka 
H-like line are about 5 eV at temperatures ~ 10 keV, 
and the natural width of the line is about 1 eV. The 
Fe lines can, however, be broadened by the velocity 
dispersion of the accretion flow in the emission region. 
The infalling material will hit the accretion shock with 
speeds ~ 6000 km s^^; if the shock is adiabatic, the 
post-shock speed will be ^ 4 times slower. Thus, the 
velocity dispersion of emission lines in an orbital-phase 
averaged spectrum could reach up to 1000 kms~^, cor- 
responding to a broadening of about 20 eV. 

Distortion of the line profiles can also be caused 
by electron (Compton) scattering. A line photon can 
gain energy through collisions with hot electrons (i.e., 
fcTc > Ec, where To is the electron temperature and 
Ec the line centre ener gy), and the fractional energy 
chan ge per scattering is (|Pozdnvakov. Sobol. fc SunvaevI 
I1977II 

For cold electrons (feTc <^ Ec), the photon loses energy 
due to recoils. For a He-like Fe line, the energy shift 
per scattering is 

Although the energy shift per scattering is small, the 
Fe line can be broadened significantly after multiple 
scatterings. 



Figure 1: Schematic illustration of the geometry 
used to model accretion onto a white dwarf, line 
emission by shock-heated electrons, and Compton 
scattering by both hot and cold electrons in the 
accretion column. 



For an mCV with a specific accretion rate of ~ 
1 — 10gcm~^s~^, the thickness of the shock-heated 
region is Xb ~ 10^ cm, and the electron number den- 
sity is rio ~ 10^^ cm~^. The corresponding Thompson 
optical depth is r ~ 0.1. The mean number of scat- 
terings is Ai" ~ max(r^,r) (Rybicki & Lightman 1979). 
Thus, one in every ten line photons would encounter 
an electron before leaving the post-shock region. The 
effective Doppler broadening due to Compton scatter- 
ing by hot electrons in this downstream region is about 
5% of the line-centre energy. 

Compton scattering by the cold, upstream elec- 
trons in the pre-shock region can, on the other hand, 
cause much more substantial line broadening. Despite 
the fact that the electron number density in the up- 
stream fiow could be an order of magnitude lower than 
that in the downstream, post-shock region, the lin- 
ear size of the accretion column between the accretion 
shock and the magnetospheric coupling region could be 
more than two orders of magnitude larger (~ 10^ cm) 
than the thickness of the shock-heated region. The ef- 
fective line-of-sight electron scattering optical depth r 
is therefore anisotropic, with values larger than unity 
at certain viewing angles. We would therefore ex- 
pect observable Compton recoil signatures in the X-ray 
emission line spectra of mCVs. 

Note that line broadening due to resonance scatter- 
ing is a factor of (me/mi)^''^ smaller than broadening 
due to Compton scattering. For the Fe H- and He-like 
lines, this factor is ~ 1/320 (Pozdnyakov et al. 1977). 
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3 Modelling Compton Scatter- 
ing in the Accretion Column 

3.1 Line emission and scattering re- 
gions 

We investigate Compton scattering effects on Fe Ka 
lines in mCV accretion columns using Monte Carlo 
simulations. We use line centre energies of 6.4 keV, 
6.675 keV, and 6.97 keV, corresponding to transitions 
of neutral, He-like, and H-like Fe, respectively. The 
neutral and Lyman-a transitions are actually doublets, 
with energies 6.391/6.404 keV and 6.952/6.973 keV 
and branching ratio 1:2, and the He-like transition has 
both resonant and intercombination components, and 
dielectronic satellite lines are also associated with the 
H-like and He-like lines. These subcomponents remain 
unresolvable with current X-ray instruments, and in 
this work, we simply use the approximate line cen- 
tre energies in order to demonstrate the importance of 
Compton scattering effects on the line profiles. 

In the simulations, the line photons are emitted 
and scattered in the accretion column, which is as- 
sumed to be a cylinder with a circular cross section. 
The cyclinder is divided into a geometrically-thin, down- 
stream shock-heated region, and a semi-infinite, up- 
stream pre-shock region. The radius of the accretion 
cylinder is assumed to be 0.1 -Rw. The line photons 
are emitted in the shock-heated region, and scattering 
occurs in both the post- and pre-shock regions. The 
geometry of the model is shown in Figure 1. We note 
that the X-ray emission from the post-shock accretion 
column illuminates the white dwarf atmosphere, which 
could cause Fe Ka fluorescence and additional Comp- 
ton recoil effects. Such effects can be modelled using 
a Monte-Carlo technique. In this study, however, we 
focus only on the scattering of the lines emitted from 
the shock-heated plasma. Our future study will ad- 
dress this issue properly. 

The post-shock accretion fiow is cooled by emitting 
bremsstrahlung X-rays and cyclotron optical/infra-red 
radiation. Line cooling and Compton cooling are gen- 
erally unimportant in most of the post-shock region 
and are thus neglected in this study. The velocity, 
density and temperature profiles in the shock-heated 
region are assumed to be nonuniform, in contrast to 
previous studies, and are calculated using the hydrody- 
namic model (with bremsstrahlung and cyclotron cool- 
ing) described in Wu, Chanmugam & Shaviv (1994). 
Although resonant scattering in the post-shock region 
can be important, we neglect this process in the simple 
model presented here, although resonance line trap- 
ping is expected to enhance the overall effects of Comp- 
ton scattering. 

We note that the ionisation structure and the line 
emissivi ties in the accretion column can be calculated 
directly tWu. Cropper, fc Ramsavl200lf) . In this study, 
however, we simply assume that the H-like and He- 
like lines are emitted from a particular stratum in the 
column and the 6.4 keV Fe line arises from the bot- 
tom stratum just above the white-dwarf surface. A 
more detailed study of scattering effects with a more 



self-consistent model of the spatial variation in line 
emissivities and line-of-sight viewing angles will be pre- 
sented elsewhere (Kuncic et al., in preparation). The 
upstream pre-shock region is assumed to be filled with 
cold (fcTo <S 1 keV) plasma accreting at the local 
free-fall speed vs = [GM^/{R„ + z)]^^'^, with den- 
sities determined by the accretion rate m, viz. — 
rn/{ij,muvs). 

The parameters that determine the structure of the 
post-shock region are therefore the white-dwarf mass 
Mw, the white-dwarf radius 7?w, the specific accretion 
rate m, and the ratio of the efficiency of cyclotron cool- 
ing to bremsstrahlung cooling es (see Wu et al. 1994; 
Saxton, Wu, & Pongracic 1997). While the white- 
dwarf mass is a free input parameter, the white dwarf 
radius is obtained from the mass-radius relation given 
in Koester (1987). 



3.2 Monte Carlo algorithm 

Line photons are injected in a stratum, and the prop- 
agation vector of the photons are in random directions 
in the rest frame of the flow. For each photon, the ini- 
tial propagation vector is Lorentz transformed out of 
the rest frame of the emitter into the rest frame of the 
white dwarf. We then determine a (tentative) scatter- 
ing point, the distance to which is determined using a 
non-linear rejection transport technique introduced by 
Stern et al. (1995), based on the concept of a virtual 
cross-section (see Nelson, Hirayama, & Rogers 1985; 
Kawrakow & Rogers, 2001). We note that in Stern et 
al. (1995), a constant density plasma is considered. We 
construct a more general algorithm which integrates 
the mean-free-path over the spatially varying electron 
density (CuUen 2001a, b). The Klein-Nishina formula 
is used in determing the scattering cross-section. At 
this scattering point, the momentum of the electron 
is drawn from an isotropic Maxwellian distribution of 
the local temperature. The momentum vector is then 
Lorentz transformed out of the rest frame of the local 
bulk flow to the rest frame of the white dwarf. A re- 
jection algorithm is then used to decide whether to ac- 
cept or reject the scattering event, with an acceptance 
probability o"KNKei/2crTC (CuUen 2001b), where ctkn 
is the Klein-Nishina cross-section, (tt is the Thomson 
cross-section, c is the speed of light, and 14ci is the 
velocity of the electron relative to the photon. For an 
accepted event, the energy and momentum changes are 
computed using the scattering algorit hm described in 
iPozdnvakov. Sobol. Snnvaetl (TS . 

In each simulation, the photons are tagged and fol- 
lowed until they escape from the accretion column. If 
a photon impacts with the white-dwarf surface after 
leaving the accretion column, it is treated as absorbed 
and is discarded. The remaining escape photons are 
binned to construct a spectrum. The most recent ver- 
sion of the Mersenne Twister random number gener- 
ator (Matsumoto & Nishimura 1998) was used in our 
simulations. A full description of the numerical algo- 
rithm can be found in CuUen (2001a). 
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Figure 2: Simulated Fe Ka emission lines scattered by electrons in the accretion column of a white dwarf 
of mass 0.5 Mq (left column) and 1.0 Mq (right column), and with a mass accretion rate of 1 gcm~^s~^ 
and 10gcm~^s~^. Cyclotron cooling is set to zero. The solid and dotted lines correspond to cases for 
photon injection sites near the white-dwarf surface and at the accretion shock, respectively. The profiles 
in (e) and (f) are a close-up of the 6 — 7keV region, showing the 6.4 keV line emitted near the surface, 
and the 6.675 and 6.97 keV lines emitted at the shock, under the same conditions as those for (c) and (d). 
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4 Results and Discussion 



We simulated electron scattering off Fe Ka photons 
with emission line central energies 6.4 keV, 6.675 keV, 
and 6.97 keV using a total of 10* photons. We consid- 
ered two white dwarf masses and radii: Mw = 0.5Mq, 

iJw = 9.2 X 10** cm, and Mw = I.OMq, i?w = 5.5 x 
10* cm. For each mass-radius, wc considered two dif- 
ferent mass accretion rates; m — lgcm~^s~^ and 
m = 10gcm~^s~^. The 6.4 keV photons were in- 
jected near the white dwarf surface, where the bulk 
speed of the accreting material goes to zero. The 6.675 
and 6.97 keV photons were each injected at either the 
white dwarf surface or the shock, where the bulk speed 
is equal to vs/4. 

Figure 2 shows the results for when cyclotron cool- 
ing is negligible and bremsstrahlung cooling dominates 
(es = 0). Figure 3 shows the results for when cy- 
clotron emission is the dominant cooling process (with 
Ss = 100). (A more detailed investigation into the 
line emissivity profile in the accretion column, taking 
into account coUisional ionisation, is deferred for future 
work.) All the line profiles in Figure 2 show broaden- 
ing effects due to clecton (Compton) scattering. The 
broadening near the base of the profiles is due to mul- 
tiple scatterings of photons by the hot electrons in the 
post-shock region. This is most evident in all cases 
where photons are injected at the base of the accret- 
ing column (denoted by solid line profiles in Fig. 2). 
As line photons emitted near the white-dwarf surface 
have further to propagate through the accreting ma- 
terial before escaping the column, the probability of 
scattering is larger than that for photons emitted fur- 
ther above the white dwarf surface, near the shock. 
Line photons emitted near the shock (with profiles de- 
noted by dotted lines in Fig. 2) encounter fewer elec- 
trons in the post-shock region and thus, broadening 
near the base of the line profiles is weaker. The broad 
bump caused by multiple scatterings can, however, be 
washed out by the continuum emission (which is emit- 
ted mostly at the dense base, see Wu 2000; Cropper 
et al. 2000), which also undergoes the same degree of 
scattering. 

Distortions due to Compton recoil by cold electrons 
can be seen more clearly in the line profiles of photons 
emitted at the shock. The signature of Compton recoil, 
a small shoulder red-ward of the line, can be seen espe- 
cially in the high-rh cases. Figs. 2(c) and (d). This is 
because the electron number density increases linearly 
with m, so the optical depth, and hence, the scattering 
probability, increases. A comparison of Fig. 2(c) and 
(d) also indicates that the shoulder-like Compton recoil 
signature is more conspicuous in high-mass mCVs, if 
other parameters are kept constant in the model. The 
effects of Compton recoil scattering by cold pre-shock 
electrons are less likely to be washed out by contin- 
uum emission provided that the seed line photons are 
injected near the shock (not the dense base of the post- 
shock region), and that the specific accretion rate m 
is sufficiently high, as indicated in Figs. 2(d) and (f). 

Our simulations show that the down-scattered pho- 
tons of the H- and He-like Fe lines emitted near the 
accretion shock can contaminate the neutral Fe Ka 
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Figure 3: Simulated Fe Ka emission lines scattered 
by electrons in the accretion column of a white 
dwarf of mass 1.0 Mq, and mass accretion rate 
10gcm~^s~"'^, with cooling in the post-shock re- 
gion dominated by cyclotron emission. The solid 
and dotted lines are for photon injection sites near 
the white dwarf surface and at the accretion shock, 
respectively. 
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line emitted near the white-dwarf surface. This can 
be seen in Figs. 2(c) and (f), which show close-ups of 
the line profiles over the energy range 6—7 keV, with 
the 6.4 keV line emitted at the white-dwaxf surface, 
and the 6.675 and 6.97 keV lines both emitted at the 
shock. The white-dwarf mass-radius and accretion rate 
are the same as in Figures 2(c) and (d), as indicated. 

Figure 3 shows simulated line profiles with the same 
input parameters as those for Fig. 2(b) and (c), but 
with the ratio of cyclotron to bremmstrahlung cool- 
ing (£s) set to 100. Qualitatively, we expect additional 
cooling of the electrons in the post-shock region to sub- 
stantially change the thickness and the density struc- 
ture of the post-shock region and hence, the scatter- 
ing optical depth. Moreover, it also alters the spatial 
Fe Ka line emissivity profiles in the region. Our simu- 
lations in Figure 3 indicate that this enhances the over- 
all Compton scattering effects, particularly in high rh 
mCVs. Note, however, that the enhanced broadening 
near the base of the line profiles makes the recoil signa- 
tures less discernable. On the other hand, this also in- 
dicates that the Compton broadening of the base may 
be less likely than in the case of no cyclotron cooling 
to be washed out by continuum emissions and thus, 
is potentially more observable. There is some tenta- 
tive evidence for this effect in Chandra High-Energy 
Transmission Grating (HETG) spectrum of the polar 
AM Her, particularly when compared to the HETG 
spectruum of AO Psc, which is an IP and thus, has 
less cyclotron cooling (Hellier & Mukai, 2004). How- 
ever, higher signal-to-noisc observations are required 
to confirm this Compton scattering effect. 

We note that in some mCV X-ray observations, 
the Fe lines appear to have centre energies slightly 
lower than expected (Ishida 1991, see also Hellier et 
al., 1998; Hellier & Mukai, 2004). In a typical obser- 
vation a small red shoulder caused by recoil scattering 
of line photons would not be easily resolved by instru- 
ments with insufficient spectral resolution. The result- 
ing lines in the observed spectrum would thus appear 
to be slightly red-shifted. Thus, recoil scattering is a 
plausible interpretation of the observed red-shifts of 
Fe lines in mCV X-ray spectra. Another possible in- 
terpretation, suggested by Hellier & Mukai (2004) is 
that satellite lines associated with the H-like and He- 
like transitions are responsible for the apparent line 
centre shifts. These interpretations will be verifiable 
only when spectra with sufficiently high spectral reso- 
lution and signal-to-noise ratio to resolve the fine sub- 
structures in the line profiles have been obtained (e.g. 
by Astro-E2 and other future X-ray satellites such as 
Xeus and Con-X). 

5 Summary and Conclusions 



shoulder-like feature red-ward of the line. This ef- 
fect is prominent for photons emitted near the shock, 
but is not obvious for photons emitted from the dense 
base near the white-dwaxf surface. Moreover, the ef- 
fect is most noticeable for systems with a high-mass 
(~ IMq) white dwarf, and with a high mass accre- 
tion rate (~ lOgcnr^s^^). (2) The 6.4 keV neutral 
Fe line can be contaminated by the photons down- 
scattered from the higher energy H- and He-like lines. 
(3) Lines emitted from the dense base of the accretion 
column experience multiple scatterings, and a very 
broad bump is formed at the base of the line profiles. 
However, the broad bump is unlikely to be distinguish- 
able from the contirmum emission, which is predomi- 
nantly emitted from the dense base of the column and 
which consequently suffers a similar degree of scatter- 
ing. The possible exception to this is when cyclotron 
cooling is sufficiently strong to enhance the density 
structure in the accretion column, thereby enhancing 
the effects of Compton scattering in the post-shock re- 
gion, and the prominence of the resulting broadening 
at the base of the line profiles. 
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